Modulated electronic configurations in selectively doped multilayered nanostructures 
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A simple theoretical model is proposed to describe the recent experimental results on formation 
of induced superconducting state and anomalous tunneling characteristics in selectively doped mul- 
tilayered nanostructures based on La2Cu04 perovskite. In particular, it is shown that the structure 
composed from the nominally non-superconducting (undoped and overdoped) layers turns to be su- 
perconducting with superconductivity confined to narrow regions near the interfaces, in agreement 
with the experimental observations. 
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Recent experiments by Bozovic et al [1] provided a 
new intriguing insight on electronic properties of nanos- 
tructured perovskite systems. Using thorough epitaxy 
techniques available in Brookhaven National Laboratory 
[2-4], they were able to selectively introduce a well con- 
trolled level (including zero) of Sr dopants into each 
particular La2Cu04 layer (along the c-axis) and then 
observed unusual electronic characteristics of the com- 
posite structures. For instance, a stack of 15 alternat- 
ing (La2- a; Sr 2 ;Cu04) 4 (La2Cu04) 2 blocks with x = 0.45, 
that is alternating overdoped [5] and undoped (both sep- 
arately non-superconducting) layers, revealed supercon- 
ductivity with the critical temperature T c = 30 K [1]. 
The authors interpreted this behavior as an evidence for 
carriers derealization beyond the nominally doped re- 
gion of the multilayered system. Below we propose a 
very simple theoretical model permitting a qualitative 
and semi-quantitave explanation of such derealization 
effect. 

The euristic basis for the model is the assumption that 
the collective electronic states in the multilayered sys- 
tem are superpositions of almost uncoupled (because of 
a very slow c-axis hopping t c ) planar states in each jth 
La2CuC>4 layer, formed by the fast afo-hopping t a b 3> t c 
in the energy band of width W — 8t a b around the rele- 
vant atomic level and shifted by a certain local electric 
potential cfj. The latter is related to the local charge 



densities pj — e (pj 



by mobile holes with density 



Pj and ionized dopants with density Xj (where e is the 
elementary charge) , accordingly to the discrete version of 
the common Poisson equation: 
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where a and c are the in-plane and c-axis lattice pa- 
rameters and e e ff is the (static) dielectric constant that 
effectively reduces the Coulomb field in the c-direction. 
This equation would be exact for potentials in a stack 
of mathematical planes, with uniform charge densities pj 
and separation c, and should model real L^-xSr^CuO^ 
layers where pj delocalized holes and Xj localized dopants 
are distributed in different atomic planes within the pe- 
riod c of jth layer. The adopted form of purely dielectric 
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FIG. 1: Schematic of nanostructured system with periodically 
introduced dopants (light grey circles) into consecutive layers 
of La2Cu04 along the c-axis. There are only 3 independent 
values of electronic density over 6 layers in a period. 



screening is justified in neglect of c-hopping processes, 
accordingly to their above mentioned weakness. We note 
that the charge densities pj naturally vanish both in uni- 



and undoped (pj 



0) 



formly doped (pj = 
systems. 

Otherwise, the hole carrier density pj is defined by the 
respective density of states (DOS) gj(s): 



pW/2-eipj 

Pj =2 9j{e)de, 



(2) 



including the spin factor 2. Thus the role of c-hopping in 
this model is reduced to establishing the common Fermi 
level £f for all the layers. Using the simplest approx- 
imation of rectangular DOS: gj{e) = 1/W within the 
bandwidth W, we arrive at the linear relation between 
Pj and ifiji 
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Then inserting Eq. 3 into Eq. 1 leads to a non-uniform 
linear system for the densities pj : 



Pj+i +Pj-i - (2 + a)p 3 



where the dimensionless value: 
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FIG. 2: Modulated electronic configuration by the shifted 
energy bands (solid rectangles) in the nanostructured system 
by Fig. 1, calculated for x = 0.45 and localization parameter 
a = 1. The dashed rectangles indicate unshifted energy bands 
for isolated doped and undoped layers, and the hatched stripe 
marks the interval of carrier densities where superconductivity 
should exist. 



is a single material parameter of the model, describing 
the localization degree of charge density fluctuations in 
the nanostructured system (less derealization for big- 
ger a). The advantage of Eq. 4 against an analogous 
system for potentials ipj is in eliminating the Fermi level 
(doping dependent) and, notably, this system assures the 
total electroneutrality condition J2i Pj — 0- The present 
approach can be seen as a more detailed alternative to 
the phenomenological Thomas-Fermi treatment [6]. 

It is elementary to resolve Eq. 4 for the densities 
through the doping levels: Pj — J2j> fjj'( a ) x 3' ■ The 
problem is resonably simplified considering it periodic, 
then the period of n layers at given a fully defines the 
coefficients fjj'(a) for 1 < j,j' < n. For the sake 
of definiteness, let us consider the sample system like 
that in the experiment, Ref. [1], with n = 6 and 
x\ = X2 — £5 = xq = x, X3 = X4 = (Fig. 1). The 
explicit solution of Eq. 4 in this case reads: 
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(a+l)(a + 3) 



(6) 



satisfying the evident electroneutrality condition P1+P2 + 
P3 = 2x. 

Using the soft X-ray resonant scattering techniques [4] 
for direct measurement of carrier densities in the experi- 
ment, Ref. [6], yielded: pf p ps 0.33, pf v w 0.24, pf v w 
0.15. A reasonable fit to this set can be achieved from 
Eq. 6 with the choice of a = 1: pf eor ps 0.315, pf cor w 
0.27, p3 hcor w 0.135, that is within the experimental error 
of ±0.03 from the measured values. 

In order to relate these carrier densities with the exper- 
imentally defined critical temperatures, we can employ 
the phenomenological bell-like dependence: 
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FIG. 3: Calculated critical temperature T c vs doping level x 
for the (La2- a; Sr a; Cu04) 4 (La2Cu04)2 system. Three num- 
bered curves refer to the respective layers from Fig. 2 and 
arrows indicate the doping levels of 0.15, 0.36 and 0.45 as in 
the experimental systems [1, 6]. 



Tph(p) = (Pmax -P){P~ Pmin)r*, 



(7) 



with p mia = 0.07, p max = 0.2, and T* = 9000 K (this 
curve being slightly below of the commonly reported 
T c (p) in bulk LSCO [5]). Using p = pf 1001 ' in Eq. 7 
yields the value of T c w 38 K, just that observed in Ref. 
[6]. This confirms the conclusion that the SC state in 
this system is limited to the nominally undoped layers 
3,4 as represented schematically in Fig. 2. 

One can compare the fitted value of a = 1 with the 
theoretical expression, Eq. 5, using the standard values 
a w 0.38 nm, c = 1.3 nm and W w 2 eV. This suggests 
as high value of e e ff as ~ 150, however it does not seem 
unrealistic if the static c-axis polarizability for La2Cu04 
[7, 8] is enhanced by a contribution from doped mobile 
carriers. 

The situation can be further traced at varying the 
doping level x (with a supposedly constant). Thus, for 
x = 0.45 we obtain respectively: pf 1 ™ 1 « 0.395, p| heor w 
0.34, pf cor fa 0.165, and then using this ^ hcor in Eq. 7 
results in T c ps 30 K, again in agreement with the mea- 
sured value [1]. 

At least, for the nominally optimum doping level x = 
0.15, we have: p 1 = 0.132, p 2 = 0.113 and p 3 = 0.055, 
and the SC state with almost maximum T c should persist 
only in the doped 1,2,5,6 layers separated by the insulat- 
ing 3,4 layers. This agrees with the observation of blocked 
tunneling through the undoped La2Cu04 layer sand- 
wiched between optimally doped La2-xSr x Cu04 elec- 
trodes [1]. 

Furthermore, combining the results, Eq. 6, and the 
phenomenological dependence, Eq. 7, one can easily 
build a model dependence for critical temperature of SC 
transition in the given La2- £C Sr a; Cu04 — La2CuC>4 sys- 
tem vs the doping level x. As seen from Fig. 3, this 
dependence chosen as the maximum value from three 
bell- like curves: T c (p) = maxj T p h(pj(x)), has generally a 
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with Sj — Sq exp(-Kj) where 
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FIG. 4: Modulated electronic configurations for a) single 
doped layer between undoped semispaces and b) interface be- 
tween two semispaces with different doping levels. 



non-monotonous behavior with the broadest region con- 
tributed by the 3,4 layers. It should be noted that the 
SC state realized in this region may be of special interest 
since much longer lifetimes of charge carriers in the nomi- 
nally undoped layers, in similarity with the well explored 
physics of 2D electron gas (2DEG) inverse layers in semi- 
conducting heterojunctions [9]. This can be an important 
property for envisaged superconducting devices in nano- 
taylored heterosystems [10] or excitonic superconductors 
[11]. 

The model, Eqs. 1-5 can be easily extended to other 
characteristic nanostructures. Thus, inclusion into an 
infinite stack of layers with some uniform doping level x 
of a single layer with different level x + Ax will produce 
a symmetric distribution of carrier densities pj = p-j 
that obviously tend to the asymptotic value: pj^oo — > x. 
Then, considering the reduced densities Sj = pj — x, we 
obtain from Eq. 4 an infinite set of linear equations: 



(2 + a)5 - 2Si = aAx, 
(2 + a)Sj = S j+l + 



i>i, 



(8) 



with the electroneutrality condition <5o + 2X)j>i $j = Ax. 
It can be easily checked that the system, Eq. 8, is solved 



and with the most inte resting cen tral value given by 
So = AxtanliK/2 = Ax^Ja/ (a + 4). From this function, 
it follows that the greatest part of added charge den- 
sity remains at the central layer, Sq > Ax/2, when the 
localization parameter a surpasses 4/3. Though being 
somewhat higher of that used in the previous analysis of 
periodically doped system, such value can be supposed to 
describe a stronger localization for the single layer dop- 
ing. Then it can support the experimental observation of 
persisting SC state in a single optimally doped layer sand- 
wiched between undoped semispaces (x — 0, Ax — x op t) 
[12], if po falls within the range [PmimPmax]- Contrary- 
wise, a single undoped layer (xq — 0) between optimally 
doped semispaces (x = x op t = —Ax) should possess a 
lower local density po = x op t — Sq, which more probably 
goes out of [pmin,Pmax] so that this layer would pertain 
insulating. 

Another exemplary case is the interface between two 
semi-infinite stacks of layers with different uniform dop- 
ing levels Xj ■ — x at j < — 1 and xj = x — Ax at j > 1, 
where the reduced densities can be defined respectively 



as Sj =x\— pj at j < 
evident symmetry Sj 
as: 



- 1 and Sj — pj — X2 at j > 1 with 
= S-j. Then Eq. 8 is reformulated 



(3 + a)Si - Si = Ax = 
(2 + a)Sj=Sj + i + Sj-i, 



given by Sj 
k as above 



xi - x 2 , 

3 > 2, (9) 

= Siexp[-n(j - 
and with 5% = 



and its solution is 
1)] with the same 

Ax/ 2 - a/2 + ^a(a + 4) 

At least, combining the previous cases can serve to 
explain the "giant proximity effect" observed in a thick 
underdoped layer sandwiched between optimally doped 
electrodes [13]. 

In conclusion, a simple electrostatic model model com- 
bined with 2D electronic band spectrum is used to semi- 
quantitatively explain the recent experimental findings in 
La2-KSr2;Cu04 multilayered systems, selectively doped 
with precision to single atomic layer. Exact solutions are 
found for local charge densities p in conducting CuC>2 
planes, for a number of periodic and non-periodic dop- 
ing configurations, permitting agreement with the exper- 
imentally defined p's and SC transition critical tempera- 
ture T c . The model can be used for effective designing of 
new SC (including Josephson tunnel) systems. 
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